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Abstract 

The evolution equation for the propagator of the quantum system in the optical probability 
representation (optical propagator) is obtained. The relations between the optical and quantum 
propagators for the Schrodinger equation and the optical propagator of an arbitrary quadratic 
system are found explicitly. 
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1 Introduction 

Since the beginning of quantum mechanics there were attempts at understanding the notion 
of quantum states in terms of a classical approach [TJ [2j [31 HI |5]. Within the framework of this 
^ ■ approach, a number of the so-called quasiprobability distribution functions, such as the Wigner 
c3 ■ function [5J, the Husimi function [7j, and the Glauber-Sudarshan function [U], later on unified 
into a one-parametric family [10], were suggested. But due to the Heisenberg principle, in contrast 
to the classical probabilities, all these quasiprobabilities do not describe distributions of measurable 
variables. 

A formulation of quantum mechanics that is very similar to the classical stochastic approach has 
been presented by Moyal [11], but the evolution equation suggested by Moyal was an equation for 
the quasiprobability distribution function (the Wigner function) and not for the probability. 

Recently in [121 H3], it was suggested to consider quantum dynamics as a classical stochastic 
process described namely by the probability distribution - so-called marginal distribution function 
(MDF) (which was discussed in a general context in [10J) associated to the position coordinate X, 
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taking values in an ensemble of reference frames in the phase space. Such a classical probability 
distribution is shown to describe completely quantum states [TH [15J [T6j . Detailed analysis of the 
tomographic probability picture of quantum mechanics was presented in [T7] . 

Within this approach, the notion of measuring a quantum state provides the usual optical tomog- 
raphy approach [HI [19], EH] and its extension called the symplectic-tomography formalism [H] IS, EE] • 
Both allow one to link explicitly the MDF and the Wigner function or the density matrix in other 
representations. In this way, starting from the evolution equation for the density matrix, the evolu- 
tion equation of a Fokker-Plank type for the symplectic marginal distribution function (symplectic 
tomogram) is obtained [121 HB]- The quantum evolution equation and the energy level equation 
have been found for the optical tomogram in an explicit form in [2TJ. Such equations allows one to 
define independently the marginal distribution function. Thus it may be assumed as a starting point 
for an alternative but equivalent formulation of quantum mechanics in what we call the probability 
representation or the classical-like description of quantum mechanics. 

In [22l [231 12H [25], the symplectic propagator for a wide class of Hamiltonians and its relation 
to the quantum propagator for the density matrix was found, which establishes an important bridge 
between the probability representation of quantum mechanics and other formulations such as the 
path-integral approach. 

In the present work, we extend the optical probability description of quantum mechanics and 
verify the formalism by applying it to a nontrivial case of physical interest, which is the forced 
parametric oscillator. From the viewpoint of general formalism, we find new equations that connect 
the optical propagator of a quantum system with its integrals of motion. 

This paper is organized as follows. 

In Sec. 2, we review the optical tomography of quantum states. In Sec. 3, we find general relations 
for the optical propagator and present the evolution equation for the optical propagator with the 
initial conditions. In Sec. 4, we obtain the optical propagator of an arbitrary quadratic system and 
check the general relations by direct substitution of found propagators for simple quantum systems. 

2 Optical Tomographic Representation of Quantum States 

In this section, we give a short review of the tomographic representation of quantum mechanics 
using the so-called optical tomogram [TBI EB ED] • For simplicity, we consider the case of one degree 
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of freedom with dimensionless variables, because the generalization of all our calculations and results 
to the case of any arbitrary number of dimensional degrees of freedom is obvious. 

If p is the density matrix of the quantum state, the optical tomogram is defined as follows: 

w(X,9) = Tr{p5(X - qcosO - psinO)} = (X,9\p\X,9), (1) 

where \X, 6) is an eigenvector of the Hermitian operator qcosd + psin# for the eigenvalue X 

(q\X,0) = (2 7 r|sin^|)- 1 / 2 exp L ^i^l \ . (2) 

I sin 9 J 

In terms of the Wigner function [6], the tomogram w(X, 9) is expressed as 

w(X,9) = — — f W(q,p)e^ x - qcosd ~ psine) dqdpdr]. (3) 
(2n) 2 J 

This relation can be reversed using the symmetry property of the optical tomogram 

w(X,9,t) = w((-l) k X,9 + 7rk,t), k = 0, ±1, ±2,--. (4) 
After some algebra, we arrive at 

7T +00 +00 

W{q,p) = —J d9 J J w(X,9)\ V \ e iv(x - qcos9 - psind) d V dX. (5) 

— oo — oo 

From (jHJ), using the relations between the Wigner function and the density matrix p(q,q') in the 
coordinate representation 

W(q, p) = J p(q + u/2, q - u/2) e - ipu du, (6) 
P(q, = fw ( L ^-,p) e ip ^dp, (7) 



2tt J \ 2 

we can write the relations between the optical tomogram and the density matrix in the coordinate 
representation as follows: 

w( x, e) = I J P (, + s - ^) ^4, d,, (8) 

p{o.i l') — J d6> J w(X, 9) \r]\ exp < irj I X — - - - cos 9 j > 5(q — </ — rj sin 0) d# cf?7 dX. 



(9) 



Thus, the tomogram w(X,9) contains all information on the quantum state. 
The quadrature statistics can be obtained from the optical tomogram 



(X n )(9) = J X n w(X,9) dX. (10) 



3 Quantum Evolution and Optical Propagator 



The quantum evolution equation for the optical tomogram has been found in explicit form in 
2T] . In the one-dimensional case, for the Hamiltonian H = (p 2 /2) + U(q) this equation reads 
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The evolution of the probability distribution w(X, 9) can be represented as the following integral 
relationship: 
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w(X,9,t) = J 69 J dX Il(X,9,t;X',6',t')w(X',6',t'), t>t', 



(12) 



where IT(X, 9, t; X' , 9', t') is the Green function of Eq. (llip . We call this function the optical propa- 
gator. In the probability terminology, it can be interpreted as the classical probability density of the 
system's transition from the initial position X' in the ensemble of reference frames of the classical 
phase space to the position X. 

Due to (H]), the optical propagator satisfies the symmetry property 



II (X, 9, t; X', 9', t') = II ((-l) fc X, 9 + nk, t; X', 9', t' 



t > t'. 



(13) 



Taking into account the physical meaning, the propagator satisfies the nonlinear integral rela- 
tionship in the case of N consecutive moments of time tk (k = 1, N) between t x = t in and tf = tjy 



. fN-l \ N-l 

n(Xjv, 9n, tf] Xi, 6*i, ti n ) = n(Xfc + i, 9k+i, tk+i] Xk, 9k, tk) > Y\_ dXkd9k- 

J { fe=l J k=2 
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If one takes in this relation tf — t in = Nr, tj, = t in + kr, in the limit case r — > 0, N — > oo, the 
expression for the optical propagator can be obtained in terms of the functional integral. 

We establish the connection between the optical propagator and the quantum propagator (Green 
function) for the density matrix p(q,q',t). 

For the pure state with the wave function \l/(g, t), we have 

p(q,q',t) = y(q,t)y*(q',t). 

Since the wave function at the time moment t is related to the wave function at the time moment t' 
by the Green function G(q, q', t, t') of the Schrodinger equation 



*( 9> t)= J G{q, q, t, t')V(q, t')dq, t > t' , 



for the density matrix we can write 

P{V, V', t) = J K(q, q, t; q, q, t')p(q, q, t') dqdq, t > t' , 

where K(q, q', t; q, q', t') = G(q, q, t, t')G*(q', <f , t, t') is the quantum propagator for the density matrix. 
Obviously, the quantum propagator satisfies the initial condition 

K(q, q', t'; q, q, t') = 8(q - q)5(q - q). (15) 

Usually the quantum propagator is taken to be zero at t < t'\ 

K(q,q',t;q,q',t') = 0, t<t'. (16) 

In view of the relation between the density matrix and the optical tomogram (Q, we obtain 

K(q, q', t; q, q\ t') = ^ J n (x, 6, t; l -{q + q') cos 6' + p' sin 9', 6', S) e' ip '^ 



xe 



ir)[X-((q+q')cos0)/2] 



(2 

\r)\6(q -q'-v sin 9) drj dp' dX d9 d9'. 



(17) 



Thus, if we know the propagator for the optical probability distribution w(X, 9), we can find the 

quantum propagator for the density matrix of the quantum state. 

Formula ( 1T71) can be reversed as follows: 

rr ,^ r n 1 f T , / fcsin# fcsin^ 

n(X, 9, t- X', 9', t') = — J K (q + 

r , 2q + f] sin 9' 



q 



,t;q! + 7]sm9',q',t' 



x exp < if] [X - 
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cos9' ] \ exp {— ik(X — qcos9)} \r)\ dkdrjdqdq. [\i 



After substituting the initial condition ( TT5|) for the quantum propagator into expression (|T8|) . we can 
obtain the initial condition for the optical propagator, namely, 



U(X, 9, X', 9', t') = 5{X cos{9 - 9') - X')5(sin(9 - 9')). 

Obviously, this condition satisfies the symmetry property (|T3|) . 

In view of ( !T2l) . taking into account the initial condition ( !T9l) . and defining 

Tl(X,9,t;X',9',t') = 0, at t < t' , 



(19) 



from ( ITT]) we can obtain that the optical propagator satisfies the equation 
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+ X cos 9 + i 



sin 9 d 
2 dX 



Il(X,9,t;X',9',t') 



5(X cos(9 - 9') - X')8{sm{9 - 9'))8{t - t'). 



(20) 



It is worth noting that relations f lT7|) -f lT8|) are the extensions of known results [26] for the sym- 
plectic tomogram 

M(X, fi, u, t) = J W(q, p, t )e lk{x - m - vp) dk dq dp, (21) 

in the case where the symplectic propagator is defined as 

M(A~, //, v,t) = J n(X, v, t; X', //, i/; t')M{X' , //, z/, t') dX dp! dv' . (22) 

The relation between the symplectic tomogram and the quantum propagator reads 

K(q,q,t;q,q,t) = —J — e*p\tY — + 

x H(Y, n,q — q ,t;Y' , //' ,u';t') d\id\i dY dY' dv 1 . (23) 



For the inverse relation we have 



niv a v' ' * +>\ 1 /' / -q-j/fy q + q' \ , . y/ . ;2g-^ 

n(A,/x, i/,t;A ,n,u;t) = -^pi^y exp j-z— — I A - /i— ^— I + - z/i 

x K(q,q' ,t;q,q — v',t') dqdq' dq. (24) 



Using the homogeneity property of the symplectic tomogram 

M (AX, A/i, Az/, t) = |A| _1 M(X, /i, z/, t), (25) 
we can rewrite this relation as follows: 

I1(X, //, i/, i; A', /i', i/'; t') = ( 2 J)2\ U \ J eXp { V ( X ' ~ 2 ~ ) } 

— i I X — /i — - — j > — u)K(q, q , t;q,q — fi,t) dq dq dqdv\. (26) 

Also the property (|25p enables us to find the relation between the symplectic and optical propagators: 
n opt (X, 9, t; A', &, t') = J n s (A, cos 9, sin G, t; rX', r cos 9', r sin 0'; t')\r\ dr. (27) 

4 Optical Propagator for Quadratic Systems 

As an example, we consider the system with quadratic Hermitian Hamiltonian 

H=^(QBQ) + CQ, 
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where Q = (p,q) is a vector operator, B is a symmetric 2x2 matrix, and C is a real 2- vector 
dependent on time. As known [271 128] , the system has linear integrals of motion: 

i(*)=A(t)Q + A(t), (28) 

where the real symplectic 2x2 matrix A(i) and the real vector A(t) satisfy the equations 

A = iABa y , A = iAa y C, 

with the initial conditions A(0) = 1 and A(0) = 0. The integrals of motion are defined by the 
following equation: 

d t I(t) + i[Hj(t)] = 0. (29) 

For the quadratic systems under consideration, any integrals of motion can be expressed as 
functions of two operators: A(t) = UaU~ l and A^(t) = Ua)U~ l , where U is the evolution operator 
and 

A{t) = -^(e(t)p - i(t)q) + P(t), A\t) = -A=(e*(t)p - e*(t)q) + (3*(t), (30) 
with e satisfying the equations 

e + co 2 (t)e = 0, ee*-e*e = 2i, (31) 

with the initial conditions e(0) = 1 and e(0) = i. The function (3(t) is defined from fT29l) as 

t 

m = -^Jdt'e(t')f(t'), (32) 



o 



and the integrals of motion ff28|) are expressed from A and A^ as follows: 



i(*)=(i P ) = f^l. i(t = 0)=( P ). (33) 



L <1 

The matrix A and the vector A read 



I ] \ A+At 



lfe + e* + * J_ ( i(/3 - f3*)\ 

2\i(e-e*) -i(£ -£*))' V2~\f3 + f3*)' 

The knowledge of integrals of motion ( 1331) allows one to find the Green function (or quantum prop- 
agator) for the Schrodinger equation of the system (see [27]). 

Now we expand the method of integrals of motion in order to find the optical propagator. Taking 
into account the initial condition (JT9~]) . we can write the system of equations (we put t' = 0) 

qU(X, 6, X\ 6', t = 0) = qU(X, 6, X'9', t = Q), p[l(X, 9, X' , 6', t = 0) = -p'U(X, 6, X'9', t = 0), 

(34) 



where p and q are the quadrature operators in the optical tomographic representation 

n -l 

n 

q = sin 9 



dX 



d9 +XC0s9+ 2 Sm9 dX> 



V 



— cos 9 



d_ 

dX 



-i 



A + Xsin^-lcos^. 



(35) 



Here, the primes mean the action to the primed variables of the propagators. 

During the time evolution (or formal action of the evolution operator to nonprimed variables), 
the system (I34p evolves to the following one: 

I q (t)U(X,9;X',9',t) = q'n(X,9;X',9',t), (36) 
l p (t)U(X,9;X',9',t) = -p'U(X,9;X',9',t), (37) 

where the operators I p and I q in the optical tomographic representation correspond to the operators 
(E3D: 

" d 



I. 



l -[cos9(e - £*) + sin^(£ - e*))-^ + % - 
-(s-s*) (cos^-sm.(l + X^ 



dX 



-1 













-sin^-cos^ 1+4 



% -[cos9{e + e*) + sin^(£ + £*)] A + I 



- e + e* 



cos^— — sin i 
39 



+ 

_d_ 

dX 

+ 



V2 ' 



- -1 














(e- 











sin o 1 — — cos i 
39 



1 + X 



dX 



V2 



One can solve the system of equations (1361) . ( 1371) by the following method. 

First, from (1361) the dependence on nonprimed variables is completely obtained, considering the 
primed variables as parameters, then from fl37j) the dependence on primed variables is found. The 
time-dependent factor is obtained by the substitution of the propagator to the dynamic equation 
( gPD with the initial condition ( EES]) . 

In our case, Eq. (120|) has the form 

il- (cos 2 9 + co 2 (t) sin 2 9)d e Tl + (1 - Lu 2 (t)) sin 9 cos 9(1 +Xd x )U + f(t) sin 00*11 

= 5(t)5(X cos(9 - 9') - X')5(sin(9 - 9')). (38) 

The described procedure enables us to find the optical propagator for our quantum system, but 
we can also find it from the results of [26] where the symplectic propagator of the parametric-driven 
quadratic system was found, 



U(X, fi, u) X', fi', u') = 5(X -X' + AfA- 1 A)5(Af' - A/"A _1 )e(*), 



(39) 



where A" and A/ 7 are vectors A" = (v, p) and A/ 7 = (i/, //), and Q(t) is a Heaviside step function 

0(f) 



1, at t> 0, 
0, at t < 0. 



With the help of (1271) . we can write 

U{X, 9- X', 9', t) = 5{Af^X + ATA" 1 A) - X'MAjj} 8{M[ - AVVA7 1 /NAj?}Q(t), (40) 

where the vectors A/" = (sin 9, cos 6 1 ) and A/ 7 = (sin cos 9'), now A7 1 and A 7/ are the first and 
second columns of the matrix A" 1 , inverse to A, and the indices 1 and 2 are related to the first 
and second coordinates of the vectors. In view of the time-dependent functions e(t) and (3(t), the 
propagator H(x, 9; x', 9', t) is expressed as follows: 

Ilpf , 9- X\ 6', t) = 8[X cos 9' - X'{sm 9(e* + e) + cos 9(e* - e))/2 
+ cos 9' (sin 9 (e* (3 + ifi) + cos9(e*(3 + ep*))/V2] 



x8 



„.sxa9(i* — e) + cos9(e* — e) . n . 
1 cos 9 — ■ — ■ — sin 9 



©(*)■ 



sm9(e* + e) + cos9(e* + e) 

Thus if the integrals of motion are known, e.g., the matrix A(£) and the vector A(£) are known, then 
according to formula f)40p . the optical propagator is known. In the particular case for a free motion 
(with the Hamiltonian H = p 2 /2), we have [271 12H] 

Po(t)=p, q (t) = q-pt, A(t) = 0, A(t) = (_ t °J 

and the propagator of free motion reads 

Ilf (X, 9; X', 9', t) = 5(X cos9' - X' cos9)5(cos9'(t + tan9) -sin9')0(t). (41) 

For the harmonic oscillator with the Hamiltonian H = p 2 /2 + q 2 /2 

... / cost sint \ 

A (t) = • 

\-sint cost/ 

from ( I40p we obtain 

n os (X, 9; X\ 9', t) = 5(X cos(9 -9' + t)- X') 5{sm(9 -9' + t))G(t). (42) 
The evolution equation, in this case, is very simple, 

(d t - d e )U(X, 9- X', 9', t) = 5{t)5{X cos{9 - 9') - X')5{sin{9 - 9')), (43) 
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or for the optical distribution 

(d t -d e )w os (X,e,t)=0. 

It is easily to check that (jUl) and f|42|) satisfy the nonlinear expression (fT4|) . and substitution of 
ET|) and fj4~2]) into f|T7|) provides known results for quantum propagators (see Appendix) 

1 



K{{q,q\q,q\t) 
K s(q,q',q,q,t) 



2nt 



cxp 



^2 



2t 
exp 



(?-?) 
i 



2t 



(q - q 



©(*); 



;(9 2 



</ 2 + <f -g /2 )cott 



sint 



(99 - 9 9 



(44) 

0(t). (45) 



, respectively. Thus, we have checked 



27r| sint| 

Substituting these expressions to f[T5j) gives us (jHJ) and 
the expressions obtained in the present paper. 

Note that during the procedure of dimension restoration, the Planck's constant % appears neither 
in the dynamic equation for the propagator and the marginal distribution of quadratic system nor in 
the general expression for the propagator; in this sense, the dynamic equation and the propagator are 
classical (if, of course, u and / do not depend on h). However, the initial condition for the dynamic 
equation of marginal distribution can contain the Planck's constant, i.e., to be quantum. 



5 Conclusions 

To summarize, we point out the main results of this study. 

We obtained the relations between the propagators for the optical tomogram and for the density 
matrix of quantum systems. We presented the evolution equation for the optical propagator along 
with the initial conditions. We showed the correspondence between the optical and symplectic 
propagators. As an example, we derived the optical propagators of arbitrary quadratic systems and 
checked our general expressions by direct substitution. The importance of the expressions written 
here for optical tomographic representation is connected with the fact that namely these tomograms 
are measured in the experiments on photon states where their characteristics are studied. The 
propagators and the evolution equation found in our work provide the possibility of monitoring the 
system quantum states in the process of time evolution. 
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6 Appendix 

Substituting fj42l to f|T7|) we can write 

/ s(xcos(9-9' + t) --(q + q')cos9' -p' sin 9') 5 (sin(9- 9' + t))exp(-ip'(q 
(2n) 2 J \ 2 / 



x exp ( ir](X — -(g + q) cos9) ) \r]\5(q — q — r]sin9)dr] dp' dX d9 d9' 



2 



(2tt 



l — |5(z(-l) fc -^(g + gO(-l) fc co S ^-^(-l) fe sin^)exp(-V(g-g')) 



x exp ( ir\(X — -(g + g') cos#)^ |^|5(g — g' — r]sin6')d?7 dp' dX d# 
^ exp(— ip'(g — g'))M<5(g — Q — r]sin9) 



x exp |i?7 f — (g + g') cos(t + #) + p' sin(t + #) — -(g + g') cos 9^j j d^ dp' d# 
After change of variables ??cos# = y, rj sin 9 = z and calculations we arrive at (145]) . 
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